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Abstract. We extend to the case of A and early-F type 
stars our study of the transport of matter and angular mo- 
mentum by wind-driven meridional circulation and shear 
turbulence. We show that our fully consistent treatment 
of the same hydrodynamical processes which can account 
for C and N anomalies in B type stars (Talon et al. 1997) 
and for the shape of the hot side of the Li dip in the open 
clusters (Talon & Charbonnel 1998) also explains LiBeB 
observations in stars with T^ff higher than 7000 K on the 
main sequence as well as in their evolved counterparts. 

Key words: Li dip; hydrodynamics; turbulence; Stars: 
interiors, rotation, abundances 



1. Lithium abundance on the hot side of the dip 

Due to its fragility to nuclear reactions, lithium is a very 
powerful tracer of the physical conditions in stellar inte- 
riors. During the last two decades, numerous abundance 
determinations of Li for both field and cluster low mass 
stars have indeed brought to light the occurrence of trans- 
port processes of chemicals in stellar radiative zones. 

One of the most striking signatures of this (these) 
mechanism(s) is the drop-off in the lithium content of 
main sequence F-stars in a range of 300 K in effective tem- 
perature centered around 6700 K, first discovered in the 
Hyades by Wallerstein et al. (1965) and later confirmed by 
Boesgaard & Tripicco (1986). This feature appears in all 
galactic clusters older than ~ 200 Myrs as well as in field 
stars (see Michaud & Charbonneau 1991 and Balachan- 
dran 1995 for references). 

Relatively few abundance determinations are available 
for main sequence cluster stars on the hot side of this so- 
called lithium dip. The scarceness of data is due to the fact 



that lithium can be observed in slow rotators only (V sin i 
lower than ~ 70 km.sec"^), while most of the stars on 
the hot side of the dip are actually fast rotators. In the 
Hyades, Praesepe and Coma^, these objects show lithium 
abundances close to the galactic value (logiV(Li)~3.31, 
with logA^(H)=12), except for a few Li deficient Am-stars 
(Boesgaard 1987, Burkhart & Coupry 1989, 1998^, 1999). 
In the field, the occurrence of Li-underabundant main se- 
quence (or slightly evolved) stars originating from the hot 
side of the dip is not rare (Balachandran 1990, Burkhart 
& Coupry 1991). For these objects, the lithium abundance 
shows indeed a dispersion of one order of magnitude below 
the galactic value. 

Among the different processes which have been pro- 
posed to account for the lithium distribution in Pop I 
low-mass stars'^, atomic diffusion is the only one which 
can be calculated from first principles without any arbi- 
trary parameter. Its importance has been spectacularly 
demonstrated in the Sun, thanks to helioseismology (Bah- 
call & Pinsonneault 1995, Christensen-Dalsgaard et al. 
1996, Richard et al. 1996). The nature and strength of 
the abundance anomalies (for lithium, but also for he- 
lium and metals) due to diffusion alone are expected to 
change in the domain going from G to F and A stars due 
to the shallowing of the superficial H-He convective zone 
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^ which all have approximately 700 Myrs 

^ In Praesepe, Am stars present the same Li abundance as F 

stars on the hot side of the lithium dip. 

Gravitational settling and radiative diffusion (Michaud 
1986, Proffitt & Michaud 1991, Richer & Michaud 1993), mass 
loss (Hobbs et al. 1989, Guzik et al. 1987, Schramm et al. 1990, 
Swenson & Faulkner 1992), meridional circulation (Charbon- 
neau & Michaud 1988, 1990, 1991), turbulent mixing induced 
by rotation (Vauclair et al. 1978, Vauclair 1988, Pinsonneault 
et al. 1989, 1990, Deliyannis & PinsonneauH 1990, Charbon- 
neau & Michaud 1991, Charbonnel et al. 1992, 1994, Charbon- 
nel & Vauclair 1992, Chaboyer et al. 1995, Talon & Charbonnel 
1998), gravity waves (Garcia Lopez & Spruit 1991, Montalban 
& Schatzman 1994). 
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when Teff increases and to the competition between radia- 
tive acceleration and gravitational settling. It has however 
been known for a long time (Michaud 1970, Michaud et al. 
1976, Vauclair et al. 1978) that diffusion theory actually 
predicts abundance anomalies much larger than the ones 
observed, revealing the occurrence of some perturbating 
macroscopic process (es). New calculations performed by 
the Montreal group taking into account turbulent diffu- 
sion yield results in better agreement with observations, 
once more illustrating the need for an improved treatment 
of hydrodynamical processes (Richer et al. 1999a). 

Regarding the F stars in the Li-dip region, carbon, 
oxygen and boron under-abundances are indeed expected 
in the case of pure diffusion (Michaud 1986, Turcotte et 
al. 1998), but are not observed in the Hyades members 
(Boesgaard 1989, Priel & Boesgaard 1990, Garcia Lopez 
et al. 1993) nor in the Be-deficient field stars (Boesgaard 
et al. 1977). In addition, in this framework Li settles and 
remains in a buffer zone below the convective envelope 
and should be drcdged-up after the turnoff; this is not ob- 
served in M67 subgiant stars which originate from the dip 
(Pilachowski et al. 1988, Balachandran 1995, Deliyannis 
et al. 1997). This implies that another process competes 
with diffusion and leads to nuclear destruction of lithium 
inside the dip. Finally, a macroscopic process is expected 
to limit both the predicted Li overabundances due to ra- 
diative acceleration in stars with Teff ~7000 K (Michaud 
1986), and the predicted strong Li underabundances due 
to settling for stars with Teff >7200 K (Richer & Michaud 
1993), which are not observed in open clusters. The over- 
abundance problem is less crucial compared to previous 
studies with the new radiative force calculations by Richer 
et al. (1999b), who follow the strong coupling that ex- 
ists between the atomic diffusion of one element and the 
abundance of another one (Richer et al. 1997). Indeed, in 
these stars the radiative acceleration on lithium is never 
higher than gravity, but of the same order of magnitude. 
However, the behavior of the more massive stars remains 
puzzling. 

Important clues for understanding the macroscopic 
processes that compete with atomic diffusion in A and 
F stars come from lithium observations in subgiant stars 
originating from the hot side of the dip. Looking at sub- 
giants rather than giants allows to track signatures of 
mechanisms acting deep in the star on the main sequence 
while permitting to avoid further complications due to 
possible lithium destruction in more advanced stages (see 
Charbonnel 1995). Observations by Alschuler (1975) of a 
few field stars crossing for the first time the Hertzsprung 
gap indicated that lithium depletion starts earlier than 
predicted by standard dilution for stars more massive than 
2 Mq. Later on, significant lithium depletion was observed 
in a non negligible number of slightly evolved field stars 
(Brown et al. 1989, Balachandran 1990). Recently, Dias et 
al. (1999) studied the behavior of lithium in an homoge- 
neous sample of field subgiants observed by Lebre et al. 



(1999) for which Hipparcos data allowed the precise de- 
termination of both the evolutionary status and the mass. 
They confirmed that stars originating from the hot side 
of the dip present a large range of lithium abundances 
which can not be explained by standard dilution alone 
and which reflect different degrees of depletion of this el- 
ement while on the main sequence, even if its signature 
does not appear at the stellar surface at the age of the 
Hyades (see Vauclair 1991, Charbonnel & Vauclair 1992). 
In this cluster, dilution is not sufficient to explain lithium 
values in evolved stars (which have masses of the order 
of 2.2 Mq), while main sequence stars present galactic 
abundance (Boesgaard et al. 1977, Duncan et al. 1998). 
This has been observed also in open clusters with turnoff 
masses higher than ~ 1.6 Mq (Gilroy 1989). 

In Talon & Charbonnel (1998, hereafter TC98), we 
assessed the role of the wind-driven meridional circulation 
and of shear turbulence as described by Zahn (1992) and 
Maeder (1995) in the transport of chemicals and angular 
momentum in low mass stars. We showed that the shape 
of the hot side of the Li dip in the Hyades was successfully 
explained within this framework which also reproduces the 
C and N anomalies in B-typc stars (Talon ct al. 1997). 

In the present paper we study the impact of rotational 
mixing in stars hotter than 7000 K, and extend our com- 
putations up to the completion of the first drcdge-up in 
order to confront our predictions with observations in sub- 
giant and giant stars. In §2 we recall the input physics of 
our models as well as the observational constraints on the 
evolution of the surface velocity. The global characteris- 
tics of the models are discussed in §3. We then present the 
theoretical predictions for the evolution of the LiBeB and 
helium abundances, and compare our results with obser- 
vations for both field and cluster evolved stars in §4. 

2. Models including rotational mixing and atomic 
diffusion 

2.1. The transport processes 

Our models (computed with the Toulouse- Geneva evo- 
lutionary code) include both element segregation and 
rotation-induced mixing, and we treat simultaneously the 
transport of matter and angular momentum. The internal 
rotation profile evolves self-consistently under the action 
of meridional circulation, which is treated as a truly ad- 
vective process, as described by Zahn (1992) and of shear 
stresses, treated as a turbulent viscosity (cf. Talon & Zahn 
1997). As mentioned in TC98, each of these two processes 
involves a free parameter of order unity. Here, we did not 
calibrate the parameters and simply used the values which 
had been used in previous calculations (TC98) namely 
Ch = I and i^th = 2/5. We recall that the first free pa- 
rameter describes the weakening effect of the horizontal 
turbulence on the vertical transport of chemical (cf. Zahn 
1992 for details) and the second free parameter describes 
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the weakening of the restoring forces of the density strat- 
ification on the development of the vertical shear. The 
vertical turbulent viscosity Vy is then 

where N is the Brunt- Vaisala frequency and K is the ther- 
mal diffusivity. We refer to TC98 (see also Talon ct al. 
1997) for more details. 

We do treat the atomic diffusion of H, He, CNO, No 
and Mg. The diffusion coefficients are computed with Pa- 
quette et al. (1986) prescription, and we include the diffu- 
sion due to gravity and to the concentration and thermal 
gradients. The case of the LiBcB elements is peculiar, since 
on the hot side of the dip they are (partly) supported by 
radiative acceleration which is of the same order of magni- 
tude as gravity (cf. Richer et al. 1999b). We thus decided 
to compute their evolution due to the macroscopic mixing 
only. 

2.2. Constraints on the evolution of the surface rotation 

The observational data indicate that, at the age of the 
Hyades, stars on the hot side of the dip have kept their 
initial rotation velocities (due to their very shallow con- 
vection envelope), while cooler stars have already been 
spun down by a magnetic torque. To take into account the 
large dispersion of the observed V sin i in the mass range 
we consider, we computed models with different initial 
velocities: 100 km.sec^^, which corresponds to the mean 
velocity of the Hyades stars with >7000 K, and 50 
and 150 km.sec'^ (see Gaige 1993 and TC98). When the 
star leaves the main sciquence, the evolution of its surface 
rotation velocity in our models is due only to the stel- 
lar expansion and momentum redistribution, i.e., we do 
not treat magnetic braking. This has no influence on the 
light element evolution on the subgiant branch since dur- 
ing this phase the surface depletion is largely dominated 
by dilution. 

2.3. Input physics and stellar masses 

The stellar models presented here include the same physics 
as in TC98. The radiative opacities from Iglesias & Rogers 
(1996) are complemented by the low-temperature opaci- 
ties Alexander & Ferguson (1995). The thermonuclear re- 
action rates are from Caughlan & Fowler (1988). 

We chose z = 0.20 as the initial metallicity of our 
models in order to compare our results to observations in 
Pop I stars in general*. The initial helium content is 0.28. 
The relative ratios for the heavy elements correspond to 
the mixture by Grevesse & Noels (1993), and the isotopic 
ratios are from Maeder (1983). 

* In TC98 the initial mixture was suited for precise compar- 
ison to the Hyades cluster. 



We performed our computations for three stellar 

masses, namely 1.5, 1.85 and 2.2 M0. At the age of the 
Hyades, the first one has an effective temperature at the 
edge of the Li-dip, while the last one corresponds to a 
slightly evolved star. 

In these models, we assumed that no mass loss takes 
place. Indeed, the evolution of the rotation velocities of 
A and F stars indicates that the matter leaving the stars 
provides no (significant) torque to their surface. The small 
amount of mass lost at the surface thus leaves the star with 
its own angular momentum. Taking this into account in 
the numerical simulation would only marginally alter the 
result and we thus neglected it. 

3. Global characteristics of the models 

3.1. HR diagram 
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Fig. 1. Evolution of the models in the HRD from the zero 
age main sequence up to the completion of the first dredge-up 
(only two masses are shown for clarity) . The solid, dashed and 
dashed-dotted lines correspond respectively to the stars with 
initial rotation velocities of 100, 50 and 150 km.sec"'^. The bold 
tracks are followed by non-rotating models 

Fig. 1 displays the evolutionary paths followed by our 
models in the HR diagram^. In Table 1, we give their main 

® Let us notice that, on the scale of our HR diagram, the 
structural changes due to rotation are not visible on the ZAMS. 
Indeed, rotating models have lower luminosities and effective 
temperatures due to the reduction of gravity close to the equa- 
tor by the centrifugal force, leading also to longer main se- 
quence lifetimes. Here however, for the 1.5 M© rotating at 150 
km.sec"^, the difference in effective temperature on the ZAMS 
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Table 1. Characteristics of the stellar models at the end of 
the main sequence 
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t 
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log L/Lq 
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dard ones. In particular, the slight changes in the helium 

profile (see Fig. 3 and §3.3) affect the evolution. Further- 
more, since mixing keeps adding fresh fuel to the core, 
the lifetime on the main sequence increases with the rota- 
tion rate. At the turnoff. the faster models are cooler and 
more luminous than the models with slow or no rotation. 
When the star leaves the main sequence, the differences 
in opacity between the models lead to evolution at higher 
luminosity on the subgiant branch for the faster rotators 
(see §3.3). 

3.2. Transport coefficients inside the stars 




-1 -2 -30 
log(l-m/M) 



- 1 —2 —3 
log(l-m/M) 



Fig. 3. Diffusion coefficient profiles at three ages the 1.5 Mq 
model for two different velocities. The line symbols are the 
same as in Fig 2 



Fig. 2. Diffusion coefficient profiles at the age of the Hyades 
for the 1.5 and 1.85 M© models and for the three different ve- 
locities. The bold lino represents the total diffusion coefficient, 
the thin full line shows the turbulent diffusion coefficient, and 
the dashed line, the effective diffusion coefficient related to the 
meridional circulation. In all cases, turbulence dominates the 
transport of the chemicals. After the main sequence, the abun- 
dance evolution is largely dominated by the dredge-up 

characteristics at the end of the main sequence. Values for 
standard models (with no rotation) are also quoted. 

The differences of global characteristic between the 
models including mixing and the standard models is strik- 
ing. Indeed, as a consequence of mixing, the abundance 
profiles inside the star are modified compared to the stan- 

in only 30 K while the lifetime is lengthened (due to structural 
effects only) by about 15 Myrs. 



Fig. 2 shows the diffusion coefficients inside the mod- 
els which arc still on the main sequence at the age of the 
Hyades. In Fig. 3 we follow the evolution of these coef- 
ficients for the 1.5 Mq star with the lowest rotation ve- 
locities. During the first part of the evolution {t = 78 
Myrs and also t = 610 Myrs for the slowest model), tur- 
bulence has not appeared yet and the full thin line rep- 
resents only the lower limit of viscosity, namely radiative 
viscosity. Then settling of helium is possible, as can be 
seen in Figs. 5, 6 and 7. During that first period, merid- 
ional circulation advects momentum towards the interior, 
building up the differential rotation. When it becomes suf- 
ficient and for the rest of the main sequence evolution, 
turbulent diffusion dominates the transport of chemicals. 

The situation is different for angular momentum 
though. Indeed, according to observations, main sequence 
stars hotter than the Li-dip are not slowed down by a 
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Fig. 4. Abundance profiles of LiBcB inside the 1.5 M© star at the age of the Hyades and at the turnoff. The profiles of helium 
(in mciss fraction) are shown at the turnoff only. The solid, dashed and dashed-dotted lines correspond respectively to the initial 
rotation velocities of 100, 50 and 150 km.sec"^ 



magnetic torque. As we discussed in TC98, since the com- 
plete equation for the transport of angular momentum ad- 
mits a stationary solution, these stars soon reach a regime 
with no net angular momentum flux. The rotation profile 
results from the equilibrium of the advection of angular 
momentum by meridional circulation and of the turbu- 
lent diffusion due to shear turbulence^ which counterbal- 
ance each other. This stresses the need for a fully consis- 
tent treatment of the hydrodynamical processes present 
in stars. 



3.3. Abundance profiles 

In Fig. 4 we show the interior abundance profiles of the 
LiBeB elements in the 1.5 model, at the age of the 
Hyades and at the turnoff for the three initial equatorial 
velocities. In the slowest model, the abundance profiles are 
relatively steep in the region of destruction of the LiBeB 



Remember we make the hypothesis that the major source 
of turbulence is shear. 



elements by nuclear reactions. More rapid rotation leads 
to the destruction of these light elements in a larger region. 

As can be seen on the profiles of helium (Fig. 4) and 
on the evolution of its surface abundance (Figs. 5, 6 and 
7), the effect of microscopic diffusion is noticeable what- 
ever the rotation velocity. The mixing tends to counteract 
it in the faster rotators, and leads to smaller variations of 
helium in the external regions. In the central regions, the 
mixing also affects the helium profile, leading to longer 
main sequence lifetime when the rotation rate is higher 
(see Table 1), and to slightly larger helium core mass. In 
addition, after the turnoff the deepening convective enve- 
lope of the faster rotators encounters a region with more 
helium than in slower models and thus, with a lower opac- 
ity. These stars consequently evolve at a higher luminosity 
on the subgiant branch. 

4. Evolution of the surface abundances and com- 
parison with observations 

The evolution of the LiBeB and of the helium abundances 
at the stellar surface from the zero age main sequence 
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Fig. 5. Evolution of the LiBeB elements and of the helium mass fraction at the surface of the 1.5 Mg star as a function of its 
effective temperature. The line symbols are as in Fig. 1. The open and black squares indicate the evolutionary point at 0.7 and 
0.8 Gyr respectively 



up to the completion of the first dredge-up can be fol- 
lowed in Figs. 5, 6 and 7. We indicate the location of the 
stars at the age of the Hyades (since the main sequence 
lifetime of the models depends on their rotation rate, we 
quote both 0.7 and 0.8 Gyr). We also show the surface 
abundance variations for models computed without any 
transport process (no rotation mixing and no atomic dif- 
fusion), in order to locate the beginning of the standard 
first dredge-up. In the models with rotation, the change 
of slope around 5700-5800 K is a signature of this event. 

4-.1. Main sequence stars 

Before discussing the LiBeB elements, one has to com- 
ment the evolution of the helium surface abundance. In all 
the rotating models, the effect of atomic diffusion is never 
completely cancelled. In the very early stages of main 
sequence evolution, helium settles (cf. Figs. 5, 6 and 7). 
Latter on, the increase of the macroscopic diffusion coeffi- 
cients (cf. Fig. 3) leads to a partial dredge-up of helium. In 
the fast rotators (V=100, 150 km.sec"^), this occurs very 
early on the main sequence, and the disappearance of the 



Hcii convection zone which is necessary to explain the Am 
phenomenon is not achieved. One may then find relatively 
few chemically peculiar stars with such a high rotation ve- 
locity, and the candidates would be very young objects. In 
the slow rotators (V=50 km.sec"^) however, the helium 
dredge-up occurs later, and the stars can appear as chem- 
ically peculiar for most of their main sequence life. These 
predictions agree nicely with the observed velocity distri- 
bution of Am and normal A stars (Abt & Morrell 1995). 

At the age of the Hyades, lithium, beryllium and boron 

arc never depleted at the surface of the rotating main se- 
quence stars on the hot side of the dip by more than about 
70, 40 and 5% respectively. The larger depletion factors 
in the 1.85 Mq model compared to the 1.5 Mq one are 
due to the slightly larger diffusion coefficient and to the 
smaller convection envelope that empties more rapidly. 
More important LiBeB destruction is obtained in the 2.2 
M0 model which is already evolved at the age of the 
Hyades. 

Our predictions (obtained both in this paper and in 
TC98) are compared to the observations in the Hyades, 
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Fig. 6. Same as Fig. 5 for the 1.85 M© star 



Coma and Praesepe in Fig. 8. Let us recall that while the 
stars with Teg < 6900 K lose angular momentum from 
their surface, the hotter ones conserve their global an- 
gular momentum during their main sequence life. In the 
domain we are interested in here, the weak differential ro- 
tation leads to diffusion coefficients that prevent any im- 
portant surface lithium depletion at the age of the Hyades, 
in agreement with the observations. 

The abundance of the LiBeB elements keeps decreas- 
ing while stars evolve. At the end of the main sequence, 
the maximum expected lithium depletion in the present 
models is of about one order of magnitude, in very good 
agreement with the observations in field main sequence 
stars originating from the hot side of the dip (Balachan- 
dran 1990, Burkhart & Coupry 1991). 

Jj..2. Evolved stars 

Even if no important lithium depletion is expected at the 
surface of the main sequence stars on the hot side of the 
dip at the age of the Hyades, more lithium destruction 
occurs inside the rotating models compared to the stan- 
dard case (see §3.3). This prepares the star to the lithium 
abundance variations observed in the latter evolutionary 



phases. After the dredge-up, the LiBeB depletion is thus 
more important in the rotating models than in the mod- 
els without transport processes (see Table 2), due to the 
enlargement of the LiBeB free regions. 

The predicted evolution of the surface lithium abun- 
dance in our rotating models explains the behavior in the 
stars more massive than 1.5 Mq observed by Lcbre ct al. 
(1999) and discussed by Dias et al. (1999). Indeed, for 
these objects a lithium dispersion of about two orders of 
magnitude is seen already at Toff ~ 6300 K, i.e. before 
the beginning of the standard dilution (at Teg ~ 5600 K). 
In the same homogeneous sample, the observed dilution 
factors are higher than the ones predicted by standard di- 
lution alone. This is also seen in open cluster giants with 
turnoff masses higher than 1.5 M0 (Gilroy 1989) and in 
the large sample of field giants of Brown et al. (1989). As 
can be seen in Table 2, lithium depletion in giants down 
to factors as large as 1000 is well reproduced within our 
framework. 

The effect is especially important for Li, which is 
the more fragile element. Regarding Be and B, the post- 
dilution values decrease at most by a factor of 3 (4) and 
5 (11) respectively in the 2.2 (1.5) Mq model compared 



8 



Charbonnel & Talon - The hot side of the Li dip : LiBeB beyond the main sequence 




Fig. 7. Same as Fig. 5 for the 2.2 M© star 



to the case without mixing. From HST spectra, Duncan 
et al. (1998) determined a dilution factor of 1.0±0.3 dex 
for B in two Hyadcs giants which present a hthium dilu- 
tion factor of 2.2±0.2 (sec Table 2). These data are nicely 
explained by our rotating models. 

5. Conclusions 

In TC98, we presented stellar models including the 
most complete description currently available for rotation 
induced-mixing, taking into account simultaneously the 
transport of chemicals and the transport of angular mo- 
mentum duo to the wind-driven meridional circulation. 
Wc; showed that the shape of the hot side of the Li dip 
in open clusters is well explained within this framework, 
which also successfully reproduces the C and N anomalies 
in B type stars (Talon et al. 1997). 

In the present paper, we studied the impact of rota- 
tional mixing in A and early-F type stars on the main se- 
quence and up to the completion of the first dredge-up. We 
showed that low lithium abundances measured in evolved 
stars originating from the hot side of the dip, which can 
not be explained by standard dilution alone, can be linked 



simply to destruction in the interior of stars still on the 
main sequence. Such a destruction is not visible at the sur- 
face of the Hyades main sequence stars with Tes > 6900 
K simply because they are too young. More BeB data are 
needed to check the validity of our predictions for these 
elements. 
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Fig. 8. Lithium abundance in the Hyades (crosses), Coma 
(squares) and Praesepe (asterisks; the lithium abundances for 
the Praesepe stars with Teff ~8280 K are upper values). The 
observations arc from Bocsgaard & Tripicco (1986), Thorburn 
et al. (1993), Burkhart & Coupry (1998, 1999). Our predic- 
tions (present paper and TC98) are also shown for comparison 
(wc use logN(Li)o=3.31) : The black dots correspond to cal- 
culations performed with an initial velocity of 100 km.sec"^, 
the pluses and minuses surrounded by a circle correspond to 
models with an initial velocity of 150 and 50 km. see" ^ respec- 
tively. The models with T^s > 6900 K conserve global angular 
momentum during there main sequence lifetimes, while for the 
cooler models angular momentum is lost from the surface 
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